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Abstract. Wc develop the theory of generahzed Weierstrass a- and p-functions defined 
on a trigonal curve of genus three. In particular we give a list of the associated partial 
differential equations satisfied by the p-functions, a proof that the coefficients of the 
power series expansion of the cr-function are polynomials of moduli parameters, and the 
derivation of two addition formulae. 



Introduction 

Constructive theories of Abelian and modular functions associated with algebraic curves 
have seen an upsurge of interest in recent times. These classical functions have been of 
crucial importance in mathematics since their definition at the hands of Abel, Jacobi, 
Poincare and Riemann, but their relevance in physics and applied mathematics has greatly 
developed over the past three decades. Algebraic curves are here intended as Riemann 
surfaces, unless specified to be singular. 

The study of the simplest hyperelliptic curves, namely curves of genus two, goes back 
to the beginning of the 20th century, and these are treated in much detail in advanced 
textbooks, see for example Baker (1907) [1] and more recently Cassels and Flynn (1996) 
|17j . Not so much is known about the simplest trigonal curves, which have genus three. 
The study of modular functions of these curves was originated by Picard, and reprised 
recently by Shiga [27] and his school. In this paper we study Abelian functions associated 
with the simplest general type of curve, the general (3,4) curve. This is an (n, m)-curve 
in the sense of Burchnall-Chaundy [7] 

Our work is based on the realization of Abelian functions as logarithmic derivatives 
of the multi-dimensional a-function. This approach is due to Weierstrass and Klein and 
was developed by Baker [T]; for recent developments of the theory of mult i- dimensional 
(T-functions, see Grant [12], Buchstaber, Enolskii, and Leykin [S], Buchstaber and Leykin 
jlOj . [TT] . Eilbeck, Enolskii and Previato [TB], and Baldwin and Gibbons [S], among others. 

We shall adopt as a template the Weierstrass theory of elliptic functions, trying to 
extend as far as possible these results to the case of the trigonal genus-three curve. Let 
(y{u) and p(m) be the standard functions in Weierstrass elliptic function theory. They 
satisfy the well-known formulae 

(0.1) p(m) = -— loga(u), (p')' = 4p=^ - - ^3, p" = 6p2-i^2 
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and the addition formula, which is a basic formula of the theory 

a(u + v)a(u — v) 



We present here two addition formulae (Theorems 18.11 and 19. ip . The first of these is 
for the general trigonal curve of degree four, whereas the second is restricted to a "purely 
trigonal" curve of degree four (see (12. ip ). The first main Theorem 18.11 is the natural 
generalization of (10. 2p . The authors realized the existence of the formula of the second 
main Theorem 19.11 from [SB] . However we were not able to use that paper to establish 
our result, instead working from results by Cho and Nakayashiki [12], Grant's paper [TU] . 
p. 100, (1.6), or a calculation using [9]. The crucial part is to identify the coefficients of 
the right hand sides of these two formulae. To calculate these, we used a power-series 
expansion of the cr-function, stimulated by the works of Buchstaber and Leykin [10] for 
hyperelliptic case and of Baldwin and Gibbons [5] for a purely trigonal curve of genus 
four. 

The cT-functional realization of Abelian functions of trigonal curve of arbitrary genus 
g was previously developed in [H] and [T^]. Using these results in the case oi g = 3 we 
present explicit formulae for 6 canonical meromorphic differentials and the symmetric bi- 
differential which allow us to derive a complete set of relations for trigonal p-functions, 
generalizing the above relations for the Weierstrass p-function. 

We note that we have recently developed a parallel, but more limited theory, for purely 
trigonal curves of genus four in [6], a paper which draws heavily on the results presented 
here. It is perhaps useful to compare and contrast these two cases. As demonstrated 
in Schilling's generalization of the Neumann system [28], there are basically two cases of 
trigonal cyclic covers, the order of a related linear differential operator that commutes 
with the given one of order three being congruent to 1 or 2 modulo 3, respectively. In 
each case, the action variables of the integrable system parametrize a family of curves of 
the same type, thus the family of curves in the (3, 4)-case cannot be obtained as a limit of 
that in the (3, 5) they have different dimensions. In the present paper, we develop 

the method and prove the addition formulae, together with the characterising differential 
equations, for the former case, in that the highest power of x appearing in the equation 
of the curve is 4 (= 1 modulo 3); this corresponds to the 'base' case of the Boussinesq 
equation, the smallest-genus spectral curve of an algebro-geometric third-order operator. 
In [B] , the case where the highest power of x appearing in the equation of the curve is 5 
(= 2 modulo 3) is addressed. The differences in the two cases manifest themselves in a 
number of ways, for example the parity of the a function is different in the two cases, and 
the two-term addition formulae are antisymmetric in the genus 3 case and symmetric in 
the genus 4 case. Also the results are given for the general (3, 4)-curve here, whereas only 
for the purely trigonal (3, 5)-case in [6]. It may be possible with some work to relate the 
(3,5)-case to the (3,4)-case, but this would not be straightforward and we have not yet 
attempted this. 
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Our study is far from complete, and a number of questions still remain. One of the 
first problems still to be considered should be the explicit recursive construction of the 
(T-series generalizing the one given by Weierstrass; for a hyperelliptic curve of genus two, 
this result was found by Buchstaber and Leykin [10], who also devised a procedure to 
derive such recursions for the whole family of {n, m)-curves [10], [11]. Another problem is 
the deeper understanding of the algebraic structure of the addition theorems developed 
here, in order to generalize results to higher genera. As a pattern one can consider the 
addition formula of [S] for hyperelliptic cr-functions of arbitrary genera written in terms 
of certain Pfaffians. Also, the description of Jacobi and Kummer varieties as projective 
varieties, whose coordinates are given in terms of (derivatives of) trigonal p-functions, is 
far from complete. We hope the results we present to be the first steps towards a general 
theory of trigonal curves of arbitrary genus, as well as a tool in the study of projective 
varieties which are images of Jacobians. 

The paper is organized as follows. We first discuss the basic properties of the general 
(3, 4)-curve in Section [1], and define a restricted version of this curve, the "purely trigonal 
case" , in Section [2l In Section [3l we introduce the a function for the general curve, and 
in Section H] the Abelian functions pij and their derivatives. Section \5\ of the paper is 
devoted to the various differential relations satisfied by these Abelian functions, and the 
series expansion of the a function is discussed in Section [HI in which the result (Theorem 
16.11) is new, is proved quite constructively, and is the key for the rest of papers. Let 6'^' 
be the standard theta divisor, namely the image of the Abelian map of the symmetric 
square of the curve that we consider, in its Jacobian variety J. The basis of the spaces 
r{J, (9(n6[^])) of functions on J whose poles are at most of order n along 6'^' are discussed 
in Section [71 as a preliminary to the two main addition Theorems in Sections [S] and [HI 
respectively. The first addition theorem is a two-term relation for the general (3, 4)-curve, 
and the second a three-term relation for the "purely trigonal" (3,4)-curve. Appendix 
A has some formulae for the fundamental bi-differential, and Appendix B has a list of 
quadratic three-index relations for the "purely trigonal" case only, as the full relations 
would require too much space. The web site [13] contains more details of the relations 
omitted through lack of space. 

While Sections 1 and 2 overlap somewhat with material in [SB], we believe that the 
results are useful to make the present paper reasonably self-contained. 
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1. Trigonal curves of genus three 

Let C be the curve defined by f{x, y) = 0, where 

fix, y) = y^+{nix + Hi)y'^ + {i^L2X^ + /igx + n^)y 

— [x^ + Hzx^ + jiQX^ + jjigx + /ii2), (/^j are constants), 

with the unique point oo at infinity. This curve is of genus 3, if it is non-singular. We 
consider the set of 1-forms 

dx xdx ydx 

(1-2) UJl = —. r, UJ2 = —. r, UJ3 



fyix,y)' fy{x,y)' fy{x,y)' 

where fy{x, y) = -^f{x, y). This is a basis of the space of differentials of the first kind on 
C . We denote the vector consisting of the forms (11. 2p by 

(1.3) LO = {uJi,UJ2,UJ3) 

We know, by the general theory, that for three variable points (xi, yi), {x2, 2/2)5 and (^3, ^3) 
on C, the sum of integrals from 00 to these three points 

U = (mi,M2,M3) 

^1.4) r{^i,yi) r{x2,y2) r(x3,yi,) 

LU + tU + to 
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fills the whole space C^. We denote the points in C'^ by u and v etc., and their natural 
coordinates in by the subscripts (ui, U2, u^), {vi, V2, v^,). We denote the lattice generated 
by the integrals of the basis (11. 2p along any closed paths on C by A. We denote the 
manifold C^/A, by J, the Jacobian variety over C of C. We denote by k, the natural map 
to the quotient group, 

(1.5) K : ^ CVA = J. 

A is a lattice of the space generated by the integrals ^ uj along any loop on C. We 
define for /c = 1, 2, 3, ... , the map 

i : Sym^XC) ^ J 

(1.6) / fP^ \ 

(Pi, ■ ■ ■ , Pfc) ^ / UJ + ---+ uj] mod A, 



oo ^ oo 



and denote its image by W^'^y (W\-^^ = J for k >3 hj the Abel-Jacobi theorem.) We will 
use the same symbol u = (mi, U2, ms) for a point m G in k~ -^(W^W). Let 

(1.7) [-1](Mi, M2, Us) = (-Ml, -U2, -Us), 

and 

(1.8) := u [-i]iy[*^l 

We call this Q^'^^ the k-th. standard theta subset. In particular, ii k = 1, then (11. 6p gives 
an embedding of C: 



l:C ^ J 



(1.9) 



p 



P ^ I UJ mod A. 



oo 



We note that 

(1.10) e[2] = iy[2]^ ew^iyw, 

differing from the genus-3 hyperelliptic case in a suitable normalization [8]. If m = 
(mi,M2,M3) varies on the inverse image k~^l{C) = k,~^{W^^^) of the embedded curve, 
we can take us as a local parameter at the origin (0,0,0). Then we have (see [26], e.g.) 
Laurent expansions with respect to Us as follows: 

(1.11) Mi = |n3^ + ---, U2 = lus^ + --- 
and 

(1.12) x{u) = — + ---, y{u) = ^ + ---. 
We introduce a weight for several variables as follows: 
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Definition 1.1. We define a weiglit for constants and variables appearing in our relations 
as follows. The weights of the variables ui, U2, for every u = [ui, U2, u-^) of W^^\{k = 
1,2,...) are 5, 2, 1, respectively, and the weight of each coefficient fij in (11.11) is —j, the 
weights of X and y of each point (x, y) of C are —3 and —4, respectively. So, the weights 
of the variables are nothing but the order of zero at oo, while the weight assigned to the 
coefficients is a device to render /(x, y) homogeneous. This is the reason why fij, /xio, fiu 
are absent. 

We remark that the weights of the variables Uk are precisely the Weierstrass gap num- 
bers of the Weierstrass gap sequence at oo, whilst the weights of monomials of x{u) and 
y{u) correspond to the Weierstrass non-gap numbers in the sequence. In particular, in 
the case considered the Weierstrass gap sequence is of the form 

0, 1, 2, 37^, 5, 6, 7, 8, 9, 10,... 

where orders of existing functions of the form x^y'^, p, g G N U {0} are overlined. 

The definition above is compatible, for instance, with the Laurent expansion of x{u) 
and y{u) with respect to u^, etc. for u E W^'^\ Moreover, all the equalities in this paper 
are homogeneous with respect to this weight. 

In the next section, we use the discriminant of C. Axiomatically, the discriminant D of 
C is defined as (one of) the simplest polynomial(s) in the /i^-'s such that D = ii and only 
if C has a singular point. Here we are regarding C as a family of curves over Z. While 
no concrete expression of the discriminant is necessary for the main results in this paper, 
we put forward a conjecture based on the results of experimentation on special cases of 
C using computer algebra. 

Conjecture 1.2. Let 

Ri = rslt^ (rsltj; (/(x, y), f^{x, y)) , rsltj^ (/(x, y), /^(x, y))) , 
(1.13) R2 = rsltj^(rslt^(/(x,y),/^(x,?/)),rslt^(/(x,?/),/y(x,?/))), 

Rs = gcd{Ri,R2), 

where rslt^ represents the resultant, namely, the determinant of the Sylvester matrix with 
respect to the variable z. Then R-^ is of weight 144 and a perfect square in the ring 

^[/^l! /^45 1^2, /^5! /^8! /^3) /^6! /^9! ^12]- 

Unfortunately checking this condition directly is a computing task presenting consider- 
able difficulties due to the size of the intermediate expressions involved. We leave this as 
a conjecture and remark only that work on a full calculation is continuing. This result is 
not crucial to this paper, but we will adopt it as a working hypothesis (see Remark 16.21) . 
To continue, we define here the discriminant D of C by a square root of R^: 



(1.14) 



D = JR 
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We comment on the choice of this root in 16.21 If the conjecture is true, D is of weight 72. 
For the convenience of the reader we give -R3 

i?3 = (256/212=^ - 27/212^3' - 128/212^6' + 144/212 Ve/is' - 192/112^9/^3 + 16/ll2/i6' 

- 80/ii2/i9Ai6V3 - 4/ii2/i3V6^ + 18 flu fig 1^3^ f'G + UAflufJ^g'^ fJ^6 - 6/ii2/i9V3^ 

- Aflg'^flQ^ - 4/X9^/i3^ + flg^ fl-s"^ fle'^ + IS/ig^/ie/is - 27 flg^Y 

for the special case of /ti = /i2 = /i4 = /is = /is = (see Section [2]). 

Definition 1.3. The 2-form Q{{x,y), {z,w)) on C x C is called fundamental 2-from of 
the second kind or (fundamental second kind bi-differential) if it is symmetric, namely, 

(1.15) ^iix,y), iz,w)) = n{{z,w), (x,?/)), 

it has its only pole (of second order) along the diagonal of C x C, and in the vicinity of 
each point (x, y) it is expanded in power series as 

(1.16) fi((x, y), {z, w)) = (j^^^ + 0(1)) «' (as (x, y) {z, w)), 

where ^ and ^' are local coordinates of points (x, y) and {z,w). 
We shall look for a realization of Q{{x,y), {z,w)) in the form 

F((x, y), (2, w))dxdz 



(1.17) 

where F{{x,y), {z,w)) is a polynomial in its variables. 

Lemma 1.4. (Fundamental 2-form of the second kind) Let S((x, y), (^, ty)) be the 
meromorphic function on C x C , 



1 

;i.l8) mx,y),{z,w)) = - — ,J2y'~' 

[x - z)fy{x,y) ^ 



f{Z,W) 



W 



{Z,W)=iz,w) 



where [ ]w means removing the terms of negative powers with respect to W . Then there 
exist differentials rfj = r]j{x, y) [j = 1, 2, 3) of the second kind that have their only pole at 
00 such that the fundamental 2-form of the second kind is given a^, 

(1.19) 0((x, y), {z, w)) = (^^S((z, w), (x, y)) + ^^^^^^^j 

The set of differentials {rfi, 772, 773} is determined modulo the space spanned by the ujjS of 



^Since x and y are related, we do not use d. 
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Proof. The 2-form 

(1.20) ^J:{{x,y),{z,w))dxdz 

satisfies the condition on the poles as a function of {x, y), indeed one can check that fll.20p 
has only a second order pole at {x,y) = {z,w) whenever {z,w) is an ordinary point or 
a Weierstrass point; at infinity the expansion (11.121) should be used. However, the form 
(ll.20p has unwanted poles at infinity as a form in the {z, w)-variables. To restore the 
symmetry given in (11.151) we complement (ll.20p by the second term to obtain (11.190 with 
polynomials rjjix, y) which should be found from Ol.lSp . That results in a system of linear 
equations for coefficients of ?7j(x, y) which is always solvable. As a result, the polynomials 
r]i{x,y) as well as F{{x,y), {z,w)) are obtained explicitly. □ 

Remark 1.5. The 1-form 

is the differential of the third kind, with first order poles at points (x, y) = {zi.,wi) and 
{x,y) = {z2,W2), and residues +1 and —1 correspondingly. 

Remark 1.6. The realization of the fundamental 2-form in terms of the Schottky-Klein 
prime-form and ^-functions is given in [1], no. 272, and the theory based on the 6'-functional 
representation is developed in [18]. Here we deal with an equivalent algebraic represen- 
tation of the fundamental 2-form which goes back to Klein and exhibit an algebraic 
expression for it, that is also mentioned by Fay in [18] where the prime-form was defined. 
The above derivation of the fundamental 2-form is done in [1], around pg. 194, and it was 
reconsidered in [T^] for a large family of algebraic curves. The case of a trigonal curve of 
genus four was developed in [S], pp. 3617-3618. 

It is easily seen that the rjj above is written as 

hj(x,y) , 

1.21 r],{x,y) = ^^-^dx, j = 1,2,3, 

where hj{x, y) G Q[/ii, /i2, /i4, /is, /is, /is, /ie, /^g, Aii2][a;, y], and hj is of homogeneous weight. 

The differentials rjj are defined modulo the space of holomorphic differentials with the 
same weight, but it is possible to choose the standard rjj s uniquely by requiring that for 
each j = 1, 2, 3 the polynomial hj{x,y) do not contain monomials corresponding to non- 
gaps with bigger j. Moreover there exist precisely 2g = 6 monomials defining standard 
differentials, for more details see [11], Chapter 4. 
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In particular, straightforward calculations lead to the following expressions 

h{x,y) = -x^ 
h2{x,y) = -2xy + fiix'^ , 
(1.22) hi{x, y) = -{5x^ + (/iiyU2 - Sfi3)x + /i2/i4 + /^e)?/ + yU2y^ + 3fiix^ 

- (yU2^ + 2yU3yUi - 2/i4)x^ - (yU5/i2 + yUe/^l + fi3fJ'4)x + |/ii/x(x, y) 

The orders of monomials defining standard differentials are printed in bold: 



0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,..., 

these can be written as 3z + 4j, < z < 2, < j < 1. We remark that the last two terms in 
the definition of hi{x, y) are chosen to provide the standard differentials described above. 
The polynomial F((x,?/), {z,w)) in (11.171) is of homogeneous weight (weight —8), and is 
given explicitly in Appendix A. 

2. Purely trigonal curve of degree four 
In Section [H] of this paper, we restrict ourselves to the curve 

(2.1) C : = x^ + fisx^ + /iex^ + /igx + /ii2 

specialized from (11.11) . We also restrict results given in Appendix B to this case to save 
space. This curve is called a purely trigonal curve of degree four. Equivalently we can 
represent the curve (11.11) in the form 

4 

(2.2) C: y^ = l[{x-ak), 

k=l 

and evaluate the discriminant D according to ( 11.14p as 

(2.3) ^= n 

l<j<j<4 

The curve C is smooth if and only if ^ aj for all i,j = 1, ... 4. While we assume this 
to be case, results in the singular cases are obtained by suitable limiting process. 

For the curve (12.11) . the basis (11.21) of differential forms of first kind and the function S 
in (11.181) can be written explicitly as 

, , dx xdx ydx dx 

(2.4) ui = — , UJ2 = ^3 = ^ = 

and 

(2.5) E((x, y), (z, w)) = ^ ^^J^^^y2 ' 
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respectively. The function a{u) is defined by using these. Let 

The curve C has an automorphism {x,y) ^ {x,(y), and for u = {ui,U2,u^) G k,~^l{C), 
acts by 

(2.6) 



{dui, du2, dus) 



This action naturally induces an action on k, ^Q^^\ [k = 2, 3, . . . ), implying that the set 
O''^] is stable under the action of [(^^. 

3. The it-function 
We construct here the a-function 

(3.1) a{u) = a{ui,U2,U3) 

associated with C for m G (see also [5, Chap.l). We choose closed paths 

(3.2) ai,(3, (l^z,j^3) 

on C which generate Hi{C, Z) such that their intersection numbers are ai-aj = Pi-Pj = 0, 
tti ■ Pj = Sij. 

Define the period matrices by 



(3.3) 



- / //I 













J 



/ in 

.V V J 



i,i=l,2,3 



Vj / Vj 



jj=l,2,3 



We can combine these two matrices into 
(3.4) M = 

Then M satisfies 



LU' LU" 

r]' T]" 



(3.5) 



M 



This is the generalized Legendre relation (see (1.14) on p. 11 of 0). In particular, u' ^uj" 
is a symmetric matrix. We know also that 



(3.6) Im (cij u ) is positive definite. 

By looking at (II. 2p . we see the canonical divisor class of C is given by 4oo and we are 
taking oo as the point, the Riemann constant is an element of (|Z)^ (see |22], Coroll.3.11, 
p. 166). Let 



(3.7) 



6 :-- 



6' 
5" 
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be the theta characteristic which gives the Riemann constant with respect to the base 
point oo and to the period matrix [uj' uj"]. Note that we use 5', 5" as well as n in (13.81) as 
columns, to keep the notation a bit simpler. We define 

cr(n) = cr(u; M) = a{ui, U2, u^] M) 



(3.8) 

where 
(3.9) 



/ 1 / /—I t \ om / /—I t /—I ii\ 
cexp[—^uriuj u)i7[d\[uj u; uj uj ) 

cexp{—^uri'u'~^ *m)x 



X 

n 



exp [27rz{| '{n + 6')uj'-^uj"{n + 6') + '{n + 6')iuj''^ 'u + 5")]] , 



3 \ 1/2 



\UJ' 



with D from (11.141) . Here the choice of a root of (13.91) is explained in the remark [6l2] 
below. The series (13.81) converges because of (13.61) . 

In what follows, for a given m G C^, we denote by u' and u" the unique elements in 
such that 

(3.10) u = u'uj' + u"uj". 
Then for u, v e C^ and £ (= I'uj' + e'uj") e A, we define 

L{u,v) ■.= u{r]"v' + r]"'v"), 

(3.11) x(^) := exp[nV^{2{i'6" - £"6') + ti")] (g {1, -1}). 

In this situation the most important properties of a{u] M) are as follows: 

Lemma 3.1. The function cr(n) is an entire function. For all u G C^, £ G A and 
7 G Sp(6, Z) , we /laue 

(3.12) a{u + t,M) = x{iMu;M)expL{u + liJ), 

(3.13) a(u;7M) = a(u;M), 

(3.14) u I— (t(-u; M) has zeroes of order 1 along G'^^, 

(3.15) a{u; M) = ^ u e e^^l 

Proof. The function a is clearly entire from its definition and from the known property 
of theta series. The formula (I3.12p is a special case of the equation from [T] (p. 286 in 
the 1995 reprint, £.22). The statement (13.131) is easily shown by using the definition of 
cr{u) since 7 corresponds to changing the choice of the paths of integration given in (13.31) . 
The statements (I3.14p and (I3.15P are explained in [1], (p. 252). These facts are partially 
described also in [Hj, (p. 12, Th.1.1 and p. 15). □ 

Lemma 3.2. The function a{u) is either odd or even, i.e. 

(3.16) a{[-l]u) = -a{u) or a{[-l]u) = a{u). 
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Proof. We fix a matrix M satisfying (13.51) and fl3.6p . Tlierefore tlie bilinear form L{ , ) 
is fixed. Tfien the space of tlie solutions of (13.121) is one dimensional over C, because tlie 
Pfaffian of the Riemann form attached to L( , ) is 1 (see [23], Lemma 3.1.2 and [20], p. 93, 
Th.3.1). Hence, such non-trivial solutions automatically satisfy (13.131) and (I3.15P : while 
(I3.14P requires the constant factor to be the same, this is guaranteed by the definition 
of a and the fact that (13.90 is independent of 7. In this sense, (I3.12p characterizes the 
function a{u) up to a constant, which depends only on the fijs. Now considering the loop 
integrals for u in the reverse direction, we see that 

[-1]A = A. 

Hence u ^ a{[—l]u) satisfies (13.121) also. So there exists a constant K such that 

a{[-l]u) =Ka{u). 

Since [—1]^ is trivial, it must be = 1. □ 

Remark 3.3. In fact cr{u) is an odd function as we see in the Theorem 16. 1[ 

We need the power series expansion of a{u) with respect to ui, U2, u^. To get the 
expansion, first of all, we need to investigate Abelian functions given by logarithmic 
(higher) derivatives of a{u). We shall examine this in the next Section. 

4. Standard Abelian functions 

Definition 4.1. A meromorphic function u 1— > ^{u) on C'^ is called a standard Abelian 
function if it is holomorphic outside ^"^(O'^]) and is multi-periodic, namely, if it satisfies 

(4.1) <P(n + Lu'n + uj"m) = ^{u) 

for all integer vectors n,m and all u ^ /€~^(6[^1). 

To realize the standard Abelian functions in terms of the cr-function, we first let 

(4.2) 

* dui dvi 

for u = (mi,M2,M3) and v = (fi, f2, "ys)- This operator occurs in what is now known as 
Hirota's bilinear operator, but in fact was introduced much earlier in the PDE case by 
Baker ([3], p. 151, [3], p. 49) (see also [13] )• We define fundamental Abelian functions on 
J by 

(4-3) pij{u) = -^^AiAja{u)a{v)\^=u = -g^loga(n). 

It follows from (13.161) that these functions are even. For the benefit of the reader familiar 
with the genus one case, we should point out that the Weierstrass function p{u) described 
in eqn. (10. ip would be written as pu{u) in this notation. It is clear that they belong to 
/^(J, 26'^!). Moreover, we define 

(4-4) pijk{u) = ^Pij{u), pijM{u) = £-^pijk{u). 
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The three index p-functions are odd and four index p are even. The functions (14.31) and 
(14.41) are standard Abehan functions from Lemma IXTl Following (and generalizing) Baker 
([3], pg 151, [1], pp.49-50) (see also [9], pp.18-19, or [I2]), we define 

, , (5ijfc^(n) = -2^272 AiAjAfcA(,a-(n)a-(t;)|^=„ 

(4.5) 

= Pijkliu) - 2{pijpk£ + pikpji + piipjk)iu), 

which specializes to 



Qijkk — Pijkk — '^PijPkk — ^Pikpjk-, Qiikk — Piikk — '^Piipkk — ^p: 

Qikkk = Pikkk — Qpikpkk, Qkkkk = Pkkkk — Qp. 



2 



kk ■ 

A short calculation shows that Qijki belongs in r(J,0{2Q^'^^)), whereas pijke belongs in 
r{J, (9(40[^1)). In particular Qisss plays a key role in what follows. 

Note that although the subscripts in pijM do denote differentiation, the subscripts 
in Qijki do not denote direct differentiation, and the latter notation is introduced for 
convenience only. This is important to bear in mind when we use cross-differentiation, 
for example the pijki satisfy 

O^Pijkiiu) = -£-^Pijkni{u), 

whereas the Qijke do not. The following useful formula (14. 6 p involving fundamental 
Kleinian functions, for the case of the general curve (11.11) . was derived in [9]. It would be 
helpful for the reader to read [2], p. 377, for the case of hyperelliptic curves. The formula 
(14.61) below is proved similarly. 

Proposition 4.2. Let m G and {xi,yi), {x2,y2), {x^^y^) be Abelian preimages of u, 
i.e. 

r{xi,yi) r{x2,y2) r{x3,y3) 

U = I id + UJ + U! 

J oo J oo J oo 

with appropriate paths of the integrals. Let (x, y) be an arbitrary point on the curve C. 
Then, for each k = 1, 2, 3, the following formula holds 



(4.6) [1 X y] 



/ r{x,y) \ 




- V 


Pij ij uj-uj 




Xk 






-Vk- 



F{{x,y), {xk,yk)) 



[X - Xk) 



2 



where F(^{x, y), {z, w)) is a polynomial defined by (11.171) or ( 1A.3I1 



Proof. Using (I3.14p and relations of differntials of the second kind on C with ones of 
the third kind (see [1], p.22,£.15 and p.22,£.ll), we have an equation connecting the 
theta series appeared in (13.80 and differentials of the third kind (see [1], p. 275, £. — 11, for 
example). Then such the equation is modified to a form suitable with cr{u) and the 2- form 
fl{{x,y), {z,w)) of (I1.19p . Finally, after taking logarithm of the modified one, operating 
Q^Q^ to it gives the desired equation. □ 
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Proposition 4.3. Suppose the (xj, yi)s and u are related as in Proposition \4.'2i Let [x, y) 
be any one of [xi,yi)s. Then we have infinitely many relations, of homogeneous weight, 
linear in 

pij{u), pijk{u), pijkiiu),--- {i, j, k = l, 2,3), 

and whose coefficients are polynomials of x, y and fijs. We list the first three of them of 
lower weights as follows : 

(4.7) 

P33{u)y + p23{u)x + pisiu) = 

{p23{u) + lfilp33{u) - P333iu)) y + {p22{u) - p233{u) 
(4.8) 

+ llJ'lp23{u))x + |yUiPi3(n) + puiu) - pi33(n) = 2xy - fyUlX^ 
- + (|p33/i2 + |p3333 - |yUlp333 + 1/^1^33 + 2/iiX - §^233 + 2/i4) V 

+ (|/U2 - |/Ul^) x'^ + (-|/ilp233 + At5 + |p2333 + |p23/^2 + 1/^1^23 - |p223)a^ 

(4.9) 

+ |pl333 + |/^2pl3 + + |yLtlVl3 " |pl23 - |AilPl33 = 0. 

More equations of this type are available in [13] . 

Proof. These relations are derived from (14. 6 p by expanding (14.61) . with respect to a local 
parameter t = x~^^^, in the vicinity of the point at infinity, and comparing the principal 
parts of the poles on both sides of the relation (14. 6p . we find the solution of the Jacobi 
inversion problem. □ 

Remark 4.4. (1) In the case of trigonal curves, formula of this type was first given explicitly 
for a particular case of the curve (ll.ip in [15] . 

(2) We use in the proof of Lemma [5?T] below the first seven relations in 14. 3[ Namely, those 
of weight from —6 to —12. 

The first two of relations in 14.31 give solution of the Jacobi inversion problem (see also 

my- 

Corollary 4.5. Suppose the {xi, yi)s and u are related as in Theorem \A.'2i The solution of 
the Jacobi inversion problem is given by {xi,yi), {x2,y2) jdiT'd (^3,^/3), where these points 
are the set of zeros of the equations (14.70 . (14. 8 p for {x,y). 

We remark that the right hand side of equations (14.70 . (14. 8 p are related to the polyno- 
mials h3{z,w) and h2{z,w) defining the canonical meromorphic differentials ris{z,w) and 
r]2{z,w). Further the first equation in (4.8) is directly related to the determinant of the 
matrices constructed in [26] , using the algebraic approach developed in [21] . 
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If we take the resultant of (14.71) . (14.81) with respect to y, we find a cubic equation in x 
which can be used to substitute for in terms of lower powers of x 

= \ (3p23 + Atlp33 - Pass) + \ {P33P22 + 2pl3 + p23p333 - P33p233 " P23^) X 
+ \p33pl2 - |p33pl33 " \pl3p23 + |pl3p333- 

If we now take the resultant of (14. 7p . (14. 9 p with respect to y, we get a quartic in x which 
can be reduced to a quadratic by repeated use of (I4.10p . This quadratic in x is not further 
reducible A quadratic equation in x has at most only two solutions and u has three free 
variables. Hence each the coefficients of 1, x, x"^ of the quadratic must all be identically 
zero. Furthermore, each coefficient can be split into two parts which are even and odd 
under the reflection (II. 7p . and each of these parts must vanish. So each term of order 
higher than two in the expansion of (14. 6 P can give up to six separate equations involving 
the p functions. The simplest two arising from the resultant of (14.70 . (14.90 are 

(4.11) p222 - 2p33p233 + 2p23p333 " yU2p233 + yU3p333 + ^^lp223 = 0, 

(4.12) P23P233 - 2p33p223 + P333p22 + 2pi33 + /il(p23p333 " P33p233) = 0, 

where pij = pijiu) and pijk = pijkiu). 

5. Equations satisfied by the Abelian functions for the general 

trigonal case 

We can use the expansion of (14.60 as described in the discussion following Theorem 14.21 
to derive various equations which the Abelian functions deflned by (14.40 and (14.50 must 
satisfy. We consider flrst the 4-index equations, the generalizations of p" = — in 
the cubic (genus 1) case. 

Lemma 5.1. The A-index functions pijke associated with (12. ip satisfy the following rela- 
tions : 

P3333 = 6p33^ + /il^p33 - 3p22 + 2/iip23 - 4/i2p33 - 2/i4, 
P2333 = 6p23p33 + /il^p23 + 3/i3p33 - fJ'2p23 - 1^5 - Atlp22, 

P2233 = 4p23^ + 2p33p22 + Atl/i3p33 " /i2p22 + 2/i6 + 3/i3p23 + Atl/i2p23 + 4pi3, 

P2223 = Qp22p23 + 4/ilpl3 + /^1/Lt3p23 + l^2fJ'3p33 + 2/-i3/i4 + /i2^p23 + 4/i4p23 + 3/i3p22 
+ 2fiiflQ + /i2Ai5 - 2/i5p33, 

P2222 = 6p22^ - 2/i2/i3p23 + AtlAi2/i5 + 2/iiyU3/i4 + 24pi3p33 + 4/ii^pi3 - 4/i2pl3 " 4pi333 
+ 4/i5p23 + 2/iiV6 - 2;U2/i6 + yU3/i5 - 3/^3^33 + 12yU6p33 + 4/i4p22 
+ /U2V22 + 4/ii/i3p22, 

P1233 = 4pi3p23 + 2p33pl2 - 2/iip33pi3 - |/il^pl3 + 1/^1^1333 + |^lVl2 + 3^3pi3 
+ |yUi/i8 + |yUl/i2pl3 - /i2pl2 + /ig, 
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P1222, = 4p23pl2 + 2pi3p22 - 2^2P33P13 - 2/^8^33 " 1/^8/^2 + 1/^2^1333 + 3/^3^12 + 4/i4pl3 

+ 1/^2^13 - 2pii - V2P13 + |/il/i2pl2 + /il/i3pl3, 
P1222 = 6p22pl2 + 6yU9p33 - /i3pl333 + 4yU5pi3 + /i2^pl2 - /Lt2/i9 + 4/i4pi2 - 2/iipii 

+ 6/i3p33pi3 - 3/i2/i3pl3 + AnV3pl3 + 3/ii/i3pi2 - fJ^l 1^2 fJ'S , 
P1133 = 4pi3^ + 2p33pii - /i9p23 + 2/i6pl3 + At8p22 - /i5pl2 + |/i4pl333 + 1/^4/^8 

+ 2/i2/i8p33 - 4/i4pl3p33 + f /^2Ai4pl3 + Ail/^9p33 " fJ'lfJ'8p23 + ^l/^5pl3 

- |yUl^At4pl3 + |yWl/i4pl2, 

P1123 = 4pi2pl3 + 2p23pll + 2/i3/i4pl3 " ^3/^8^33 " 2/i5pl3p33 + /^2/i8p23 + |/i2/i5pl3 

- /i9p22 + 2/i6pl2 + |yU5pl333 + |yU5/^8 + fJ'lf^9p23 - 1/^1^5^13 + |/ilyU5pl2, 
P1122 = 4pi2^ + 2piip22 + |/il^/i6pl3 + |/^1/^6P12 + yU3/^9p33 + /^2/i9p23 + 8/ii2p33 

+ 2/i3/X4pl2 - f /i6pl333 + 4/i8pl3 - 1/^6/^8 + 4/i6p33pl3 " fJ'3fJ'8p23 + fJ'3f^5pl3 

- |/^2/^6pl3 + /i2/i8p22 + fJ'2fJ'5pl2, 

P1II3 = 6pi3pll + 6/i2Ai8pl3 - 2/i2yUi2p33 " 1^1^ I^8pl3 + 4/ii/ii2p23 + /^lyU8pl2 + /^5/^9p33 
+ At5^Pl3 - 2/i4/i9p23 + /Lti/i9pi3 - 6/i8p33pl3 " 2/^6/^8^33 + /^8pl333 " 4/^4/^12 
+ 3/X9P12 - 6/ii2p22 - /^5/^8p23 + 4/i4/X6pl3, 

P1II2 = 6pi2pll + 6/i3/ii2p33 + 3/i3/X8pl3 " 2/i6/i8p23 " /il/i8^ + 5/i2/i8pl2 + 4/i2/il2p23 

- 2/ii/ii2p22 + 4/i4/i6pl2 - f^5fJ'8p22 + /^5 Vl2 + 4/i5/ii2 - /i9pl333 - 4/ii/ii2/i4 

+ yLH^/^9pl3 + 3/ii/i9pi2 - 2//4/i9p22 + /i5/^9p23 - 4/i2/i9pl3 + 6/i9pl3p33 - 3/i8/^9, 

2 2 2 2 2 2 2 

pun = 6pii + 4/i4/i9pi2 - 8/i4 /ii2 - 2/i2 /i4/il2 - 3/i8 P22 - 2/i4/i8 + /^5 Pll - 3/i9 P33 

- 4/ii2pl333 + 24/ii2p33pl3 + 12/i5/^12p23 + l^2l^Al^bl^9 " 6/ii/i3/i4/il2 

+ /^l/^2/^5Atl2 + 2/16^8 + 2/i2V8^ " At5At6At9 - 2fir^flgpi3 + Afl^fiQpn + AfiQfIspl3 
+ 8/i2/i8pll - 6/i2/i6AH2 - 12/i2/il2pl3 + 4/iiVl2pl3 + 2/11^6/^12 + 2/i8/X5pl2 

- 6/X8/X9P23 - 12/i4/il2p22 + /i2/i5^/Lt8 + 2/ii/i4/i6/Lt9 + 1^11^51^61^8 + 12/i6/^12p33 
+ 4/ii/i9pii + 2fi3f^l4^flg + 9/i3/i5/ii2 - 2flifi3fl8^ - GfisfisfJ^g + 2fiifi2f^8f^9 

+ /i3/W4/^5At8 + 2/i2/i4/i6At8 + 2/i2/i9^- 

Proof. Many of these relations follow from the sets of equations generated from the first 
seven terms of the expansion of (14.60 as indicated in Proposition 14. 31 by a similar argument 
as that explained at the end of the previous Section. Others can be derived making use 
of derivatives of the equations in Lemma 15.51 or products of these equations with three 
index expressions pijk, working in a self-consistent way from higher to lower weights. The 
calculations are somewhat long and tedious and much facilitated by heavy use of Maple. 
Full Maple worksheets are available on request from the authors □ 
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Remark 5.2. The complete set of the four- index relations for p-functions for genus three 
was derived by Baker [5] in the hyperelliptic case only. As far as we know, the above 
relations are new, and a comparison with Baker's relations is of interest. 

Remark 5.3. With the use of (14.51) . these equations can be written in a slightly more 
compact form involving the Qijke functions. For example, the sixth equation (for P2222) 
becomes 

<52222 = -2/i2/X3p23 + /il/i2/i5 + 2/ii/i3/i4 + 4/ii ^13 " 4/i2pl3 - 4(5l333 + 4/i5p23 

+ 2fli^flQ - 2/i2yU6 + /i3yU5 " 3/13^33 + 12/i6p33 + 4/i4p22 + ("2^22 + 4/ii/i3p22- 

The importance of this switch to the Q variables is that the equations become linear in 
the Qijke and the 2- index pij. An alternative way of looking at this is that the equations 
in Lemma [5.11 have only second-order poles in a. 

Remark 5.4. The first relation in Lemma [5.11 after differentiating twice with respect to 
^3, becomes the Boussinesq equation for the function ^33 (see [HI [15]). 

Lemma 5.5. The 3 -index functions pijk associated with (12. ip satisfy a number of hi-linear 
relations {linear in both 3-index and 2 -index functions) . These have no analogue in the 
genus 1 case. For example, in decreasing weight, starting at —6 we have 

-2 P33P233 + 2 p23p333 + /^3p333 + /^2p233 - /^lp223 + P222 = 0, [-6] 
-2 P33P223 + /^1P33P233 + P23p233 " Atlp23p333 + P333p22 + 2 P133 = 0, [-7] 
-2 p23p223 + 4 P22P233 + 4/iipi33 + /i3/^2p333 + /^2^p233 

+ 4 /i4p233 - 2 /i5p333 " 2 p33p222 + yU2p222 " yUl/i2p223 " 4 pl23 = 0, [-8] 
3 Atl/i3p223 - 3 /i2/^3p233 " 24 p33pl33 + 24 pi3p333 + 12 pl22 " 12 yUipi23 

-f 12 /i2pl33 - 3 /i3p222 + 6 /i5p233 - 3 /i3^p333 + 12 /i6p333 " 6 p23p222 + 6 p22p223 = 0, [-9] 
2 P33pl23 - yUlp33pl33 + /ilpl3p333 + p23pl33 " Pl2p333 " 2 pl3p233 = 0, [-10] 
P113 + P13P223 - 2^4pl33 + P33P122 - P22P133 - P12P233 + /^8P333 - /^2pl33p33 

- /Wlpl3p233 + yU2pl3p333 + yUipi33p23 = 0, [-11] 
-Pll2 - 3 /i9p333 + P13P222 - P12P223 - 2 P22P123 - 2 /i5pl33 + fJ'lPuS + 2 p23pl22 

- At8p233 - /^2pl3p233 + 3 /i3p33pl33 - 3 /i3pl3p333 + At2p23pl33 = 0, [-12] 
8/X4P133P33 - 8/i4pl3p333 - 4/i2/i4pl33 + 2/ii/i9p333 - 2/ii/i8p233 + 2/ii/i5pl33 

-f 4/ii/i4pi23 + 4/i8/i2p333 + 3yU3pi3p233 " 3yU3p23pl33 - yUlpll2 + 3 pl2p222 
+ 4 P11P333 - 2 /i6pl33 - 3 P122P22 - 4 /i4pl22 + yU9p233 + 2 Ai8p223 " 8 p33pll3 

+ 4 P13P133 - 2 /il V113 + 2 /i2pll3 - 2 /i5pl23 = 0, [-13] 
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4 Pl23pl3 + 4 /i4p23pl33 + /^3Ai8p333 " 2 /i5p33pl33 + 2 ^5pl3p333 + /^2Ai8p233 + /i8p222 

- 4 pl2pl33 - 2 p23pll3 + 2 p33pll2 - 4yU4pl3p233 - l^lfJ'8p223 = 0, [-14] 
-/^9P222 + /ilyU9p223 + 4 pl3pl22 + 2 p23pll2 " 2 pll3p22 " 7^3/^9^333 " /^2/^9p233 

+ 2 /i5p23pl33 - 8 /il2p333 " 4 fJ^sPlSS " 4 /i6pl3p333 + 4 /i6p33pl33 " 4 pl2pl23 

- 2/i5pl3p233 = 0. [-15] 

where the number in brackets [ ] indicates the weight. 

Proof. We have already given the first two of these equations in the discussion following 
Theorem I4.2[ Some of the others follow in the same way from the expansion of (14.21) . 
Alternatively, some can be calculated directly by expressing the equations in Lemma 15.11 
in terms of pijki and pmn functions, then using cross differentiation on suitably chosen 
pairs of equations. For example the first relation above for P222 can be derived from 

d -9 _ 

7^ P3333 — -7: P2333 — 0. 

OU2 OU3 

□ 

Remark 5.6. For a fixed weight, these relations are not always unique, for example at 
weight —11 we also have the relation 

p33pl22 + 2p23pl23 + 3pll3 + /^2pl3p333-/^2p33pl33+/^8p333-2pl2p233-2/i4pl33-pl3p223 = 

These dual relations arise because in some cases the cross differentiation can be done in 
two different ways. In deriving the results in this section, it is sometimes required to make 
use of both bilinear relations at a given weight to provide enough equations to solve for 
the unknowns. A full list of the known bilinear relations is given at [13] . 

Lemma 5.7. The quadratic expressions in the 3 -index functions pijk associated with (12. ip 
down to weight —23 can be expressed in terms of {at most cubic) relations in the pmn o,nd 
P1333. For example we have the following five relations down to weight —8 : 

2 2 2 2 3 2 

P333 = P33 /^l + 2/iip23p33 + P23 + 4pi3 - 4p33p22 + 4p33 - 4^2p33 " 4/X4p33, 

P233P333 = '^fJ'3p33^ + 4p33^p23 - /^Ip33p22 - 2/i5p33 - 2/i2p33p23 + /^I^p33p23 " 2pi2 

- P22P23 + /^1P23^ + 2/iipi3, 
Pl33p333 = -|/^1P33P12 + |Ail^p33pl3 " |/^2p33pl3 + |p33pl333 - |/^8p33 + P23pl2 

+ fJ'lpl3p23 - 2pl3p22, 

P223p333 = 2/iip23p22 - 2/i2p33p22 + 2/ii/i4p23 " AtlAt5p33 + 2p33^p22 - 2/i4p22 + '^p33p23^ 
+ |/^l^Pl3 - |/^2pl3 - |/^lPl2 - 1/^8 - 2p22^ + /il/i2p33p23 + |pl333 
+ P23p33fj'3 + /^1/^3P33^ " /^2p23^ " /^5p23, 

P233^ = 4p33p23^ + 8pi3p33 + 4/i3p33p23 " 2/iip23p22 + |/ilVl3 - ^t^2pl3 + 4/i6p33 

I 2 2 4 I 2 4 4 

+ 1^1 P23 - 3/^8 + P22 - 3P1333 - 3/^1^12- 
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The expressions at lower weight quickly become very lengthy. For the purely trigonal case 
we give a list of the known quadratic expressions in the 3 -index functions up to weight —15 
in Appendix B. The full list for the general {3,A)-curve down to weight —23 is available 
at [13]. 

Proof. The relations can be found using a combination of three types of intermediate 
relations. One type is from terms in the expansion of (14. 6p . Another is to multiply one of 
the linear three-index pijk relations above by another pijk and substitute for previously 
calculated pijkpemn relations of higher weight. Yet another is to take a derivative of one of 
the bilinear three-index p^jk relations above and to substitute the known linear four-index 
pijki and previously calculated pijkpimn relations. Again, we work in a self-consistent way 
from higher to lower weights. The strategy for all the results in this section is to proceed 
down one weight at a time and to derive all the three types of relations (4-index pijke, 
bilinear 2- and 3-index, and quadratic 3-index) at a given weight before moving down 
to the next. An extra complication is that at certain weights some of the intermediate 
calculations can involve quartic terms in the pmn and P1333. It is always possible to find 
enough relations to eliminate the quartic term up to weight —23. □ 

Remark 5.8. (1) These relations are the generalizations of the familiar relation {p'Y = 
4:p^ — g2p — gs in the genus 1 theory. 

(2) For equations of weight below —23, we have not been able to find cubic expressions 
for the pijkpimn terms. We believe it should be possible to explain this using the results 
of Cho and Nakayashuiki [12], and we are currently investigating this possibility. 

(3) The calculations in this section make no use of the expansion of the a function, which 
is given in the next section. 

6. Expansion of the ct-function 

This Section is devoted to show the coefficients of the power series expansion of 
is a polynomial in /ijS. 

In the Weierstrass formulation of the theory of elliptic functions, the a-function is 
defined as the power series expansion in the Abelian variable u with coefficients depending 
on the Weierstrass parameters g2,g3, and related by certain recursive relations. The 
extension of Weierstrass theory to arbitrary algebraic curves was intensively developed 
in the 19th century and later, its development being attached to names such as Baker, 
Bolza, Brioschi, Burkhardt, Klein, and Wiltheiss. Some important modern developments 
of this theory are due to Buchstaber and Leykin [THl E] who give a construction of linear 
differential (heat-like) operators that annihilate the cx-function for any (m, ?7.)-curve. In 
the hyperelliptic case the operators are sufficient to find the recursion defining the whole 
series expansion. The exact analogue of the Weierstrass recursive series formula is known 
only for genus two, see [11], p. 68. In other cases the detailed results have not yet been 
developed, although the general method is provided in the publications mentioned above. 
Here we shall give the few first terms of the power series expansion, obtained by finding 
the coefficients of the Taylor series by using the PDEs given in Lemma 15. 1[ 
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Theorem 6.1. The function cr(u) associated with the general trigonal curve fll.ip of genus 
three has an expansion of the following form : 

(6.1) a{ui,U2,U3) = e ■ ((75(^1,^2,%) + (76(^1,^2,%) + ^7(^1,^2,%) H ), 

where e is a non-zero constant and each Cj is a polynomial composed of sums of monomials 
in Ui of odd total degree and of total weight j with polynomial coefficient in fiiS of total 
weight (5 — j). Especially, is an odd function (aee 13.21) . The first few CjS are 

C5 = Ml - Us U2^ + 2^ ^^3^ = n AH"3% - I /ilM2^, 

= 5M (/"l^ ~ 3 ^l2)uz + I /i2M3^M2^ ^'s = 3B0 (/"l^ + 9 /ig - 2 /ii/i2)M3% " 1 /is^S^Ma^, 
= 25I2O (/"l^ ^ 3 yUa) ^^^^ + (2 /i4 - /i2^ + yUl V2 + 6 /ii/is) U3^n2^ 

- ^ (4 /ii/i3 + 4 /i4 + /ig) M3M2^ + ^ /i4M3^Ml , 

'^lo = 20T60 (8AtiAt4 - 54/i2/i3 + 3/ii/i2^ + IS/iiVs + AH^ - 12 /is - 4/ii'^/i2)M3% 
+ ^ (6 yU2/i3 + 2 /iiyU4 + /ii/i2^ + /^i Va) ii3%^ 

- (4//1V3 + Ail/^2^ +4/i5 + 4/ii/i4 -2^2l^s)U2^ + ll^5U3^U2Ul, 

'^ll = ~ 6652800 (^^ /"lAi2/i3 + 27 /ii^yU2 - 72 /ig - 3 yUl*' - 24 /i2/i4 + 16 yUi ^4 - 24 /ii/ig 

+ 27 /i3' + 85 /i2^ - 4/11% - 82/^1^2')%'' + 5^ {27 fis^ + /i2^ - 6/^2/^4 

- 18//i/i2//3 + 8//1V3 - 4/ii/i5 + 6/11^4 + 12//6 + /ilV2 - ^ 1^1^ l^2^)us U2^ 

- ^ (9/i3^ - /i2^ - 4/i2/i4 -'l[lx^2[iz)uiu2^ 

+ 3^ (/il/i5 - 4/i2/i4 + /ilV4 + 3/i6)M3^Ml - \ ^tsUsU2UY, 

^12 = - 181I400 (27/ii/i3^ - 243/i2V3 - An'^ + 72/ii/i2/i4 - 31/^1^3 - 144/i2/i5 - 16/ii^/i4 
+ 6 /ii V2 - 10 /ii ^2^ + 24 /ii V5 + 4 /ii/i2^ - 72 /ii/i6 + 180 //i ^2/^3) ^i3^°^i2 
+ 21k (^^ '"3/W4 - 2 /ii//2^ + 27 //i/i3^ - 9 /i2 + /ii ^2^ + /^i Vs + 6 /ii ^2/^3 

+ 2/iiV4 + 12/ii/i6)M3%^ - ^/i3(3/il/i3 + 4/i4 + li2'')usU2 

+ llo (6 At3At4 + 2 /ii/i6 - /i2/i5 + AH V5) UsU2Ux - I (2 /il/i6 + 2 /i3/i4 + At2At5) U^U-iUx. 

Proof. We divide the proof into four parts. 

Step 1. We have aheady shown in 13.21 that all the terms are of total odd degree or even 
degree. We first show that the expansion contains a term linear in ui, so the expansion 
must be odd. 

Let B{D) be the Brill-Noether matrix for an effective divisor D of C. Then it is well 
known that (see for example [24] or [25] ) 

dimr(C, 0{D)) = degD + 1 - rankfi(D), 

where r{C,0{D)) is the space of functions on C whose divisor are larger than or equal 
to —D. Moreover, for two points Pi, P2 on C, dim r{C,0{Pi + P2)) > 1 if and only 
if the point l{Pi, P2) E O^^' is a non-singular point of 6'^^ (note that C is of genus 3). 
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By checking the Brill-Noether matrix B(Pi + P2), we see O'^' is non-singular everywhere, 
Especially ^^^(6^^^) is non-singular at the origin (0, 0, 0). On the other hand, let u and v 
be two variables on k~^(6[^1). Then we have an expansion with respect to v^: 

= a{u + v) = a3{u)v-i + \ (o-2(n) + a-iziu)) v-^ ^ , 

where cxj = da/dui, etc. Hence 

cr3(M) = a2{u) + cr33(n) = 0. 

Again by expansion 

= 0-3 (n) = o-33(0)t'3 H , 

we see that 

(X33(0) = 0. 

In summary, 

(X3(0) = ^2(0) = 0, 
so from the above arguments and (I3.14p . we must have 

ai(0) ^0. 

Hence the a-expansion must be odd. 

Step 2. Next we show that the terms of weight less than 5 vanish and C^{u) is non-trivial. 

We write all the possible odd terms up to and including terms of weight 5. Using the 
first two equations in 15.11 we can show that the coefficients of the terms of weight four 
and less are zero, and that the coefficients of weight 5 are given by those in C5 up to 
multiplication by a constant. We know from Step 1 that this constant is non-zero and we 
insert this constant into the e. 

Step 3. We now calculate the coefficents Cj, z > 5. The proof of this step is by con- 
struction (with heavy use of Maple) using the PDEs given in Lemma 15. 1[ We expand 
(7(^1,^2,^3) in a Taylor series with undetermined coefficients, keeping only odd terms. 
We do not assume that the coefficients of the expansion are polynomial in the /ij, only 
that they are independent of the Ui. We then insert the expansion into the 4- index PDEs 
for the p, and truncate to successive orders in the weights of the Uj. These give a series 
of linear equations for the coefficients, and be using a sufficient number of the PDEs we 
can always find unique solutions, as listed above. We have carried out this calculation 
down to Cig. We have omitted the details of the expressions for C13, . . . , Cig, as they are 
rather lengthy, but these are available at |13j . 

Step 4. Now consider the general term in the expansion. Set 

to be the lowest weight unknown term. Since we have already shown by construction that 
the coefficients for all weights down to —29 with respect to UjS are polynomials, we may 
assume that p + q + r ^ A. We consider the set (jj) of quadratic equations in a{u) and 



22 J. C. EILBECK, V. Z. ENOLSKI, S. MATSUTANI, Y. ONISHI, AND E. PREVIATO 

its (higher) derivatives obtained from the above, by multiplying the equations in 15.11 by 
a{uY. We take an equation 

(6.2) a{uY Qijkiiu) = ■■■ 

from (tl) such that Ui^U2^u^^ is divisible by UiUjUkUi. We have at least one such equation. 
Differentiating (16. 2p . we have an equation of the form 

/ Qp+g+r \ 

(6-3) ^{u) L ]{u) + ... = 



^duiP dui'i duf 

such that all terms are polynomial of (j{u) and its higher derivatives and such that 

Qp+q+r ^ 



[U] 



^duiP dui'i dui^\ 

is the highest derivative in (16.31) . By looking at the coefficient of the term -Ui, we have a 
linear equation of the form 

A + --- = 

over Q[/ii, ■ ■ ■ ,/ii2]. Since the other terms except A in the above equation come from 
terms of (y{u) whose weight is less than weight of ui^U2^U'i' , we see A is a polynomial in 
the /ijS by the induction hypothesis. □ 

Remark 6.2. (1) In Theorem l6.lt the constant e might be unity, another 8th root of 1, or 
some other constant. We have not been able to narrow down this result. If the case e = 1 
is true, then the determination of e reduces to the choice of roots in (I1.14p and (13.91) . The 
remaining results in this paper do not depend on this choice, or on the possibility that 

(2) The weight of ct(m) is inferred from (13.91) since the weight of \uo'\ is 5 + 2 + 1 and the 
conjectured weight of D is 72. The weight of the terms in the exponentials are all and 
the weight of c is 72/8 — (5 + 2 + l)/2 = 5 and coincides with the terms in the expansion 
of I6.1l if the weight of e is 0. 

We shall need later on the following special property of the cr-function in the purely 
trigonal case: 

Lemma 6.3. The a function associated with the purely trigonal curve (12.11) satisfies 
a{[—(]u) = —(a{u) for u G C'^ under the notation (12. ip . 

Proof. Since A is stable under the action of [Q and [—1], we can check the statement by 
Lemma 13.11 and Remark [31 □ 

7. Basis of the space r(J, C(n6[^l)) 
For notational simplicity, we denote 
(7-1) d, = £. 
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We also define 

(7.2) pf*-'' = the determinant of the (i, j)-(complementary) minor of [pijjsxs- 

We have exphcit bases of the vector spaces r'( J, 0(29'^!)) and F {J, 0(30^"^^)) as follows 
(see also [I2], Example in Section 9): 

Lemma 7.1. We have the following: 

r{j, o(2e[2])) = ci © Cpn © Cpi2 © Cp^ © Cp22 © Cp23 © Cp33 © C(5l333, 

12 © Cpiis © Cpi22 © Cpi23 
© Cpi33 © Cp222 © Cp223 © Cp233 © Cp333 © CdiQuss © C92gi333 
© C93gi333 © Cp[^^] © Cp[^2] ^ (^^[13] ^ ^^[22] ^ ^^[23] ^ ^^[33]^ 

Proof. We know the dimensions of the spaces above are 2^ = 8 and 3^ = 27, respectively 
by the Riemann-Roch theorem for Abelian varieties (see for example, [23], (pp. 150-155), 
|20j . (p. 99, Th. 4.1). Moreover, (13.141) shows that the functions in the right hand sides 
belong to the spaces of the left hand sides, respectively. For the space r{J, 0(29[^1)), pij 
and Qijki become the basis of the space from Definition 14. 11 Lemma [3?n and the arguments 
in the previous section. However these are not all linearly independent, since there are 
connecting relations, such as those given in Lemma [F!T| and the number of these relations 
is greater than the dimension of the space. Thus the problem is reduced to picking the 
linearly independent bases as a function space. It is obvious that such independence does 
not depend upon the coefficients of curve by considering these expansions around the 
origin of C^. Hence by multiplying by cr(ti)^ from the right hand side with respect to 
ui, U2, Us, and after putting all the /ij equal to zero, we see the functions of the right 
hand side are linearly independent. The authors used a computer to check this. Similarly, 
for the space r{J, (9(30^^1)), the 27 functions obtained by multiplying by a{u)^ from the 
right hand side are checked to be linearly independent by using a computer, expanding 
the given functions in the Abelian variables (cf. Theorem 16. II) to a sufficiently high power 
that independence is checked. We also see both decompositions in Lemma mi in Example 
in Section 9 of [12]. □ 

8. The first main addition theorem 

Theorem 8.1. The a-function associated with (12. ip satisfies the following addition for- 
mula on J X J : 



a{u + v)a{u — v) 



Pu{u) - pu{v) + pl2{u)p23{v) - pl2{v)p23{u) 



a{uya{vy 

+ Pl3{u)p22{v) - Pl3{v)p22{u) + | ( p33 (^i)<5l333 (f ) " p33{v)Ql333{u)) 
- {Pl2{u)p33{v) - Pl2{v)p33{u)) - /il {pi3{u) p23{v) - Pl3{v)p23{u)) 

+ I (/Ltl^ - At2) {pl3iu)p33{v) - pl3{v)p33{u)) + i/ig ip33iu) - p33{v)) 
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Proof. Firstly, we notice that the left hand side is an odd function with respect to {u, v) t— > 
([— l]n, [— and that it has poles of order 2 along (0[^1 x J)U(Jx O'^l) but nowhere else. 
Moreover it is of weight —10. Therefore, by Lemma [7.11 the left hand side is expressed 
by a finite sum of the form 

(8.2) {X,{u)Y,iv)-X,{v)Y,{u)), 

j 

where the Aj are rational functions of the /ijS with homogeneous weight, and the Xj and 
Yj are functions chosen from the right hand side of the first equality in Lemma 17. 1[ We 
claim that all the Aj are polynomial in the /ijS. Suppose all the Ajs are reduced fractional 
expressions, and at least one of the Ajs is not a polynomial. Take the least common 
multiple B of all the denominators of the AjS. Note that there is a set of special values of 
the /ijS such that B vanishes and the numerator of at least one Aj does not vanish. After 
multiplying the equation "lhs"= (18. 2p by B a^uYa^v)"^, and taking the fitS to be such a 
zero of B, we have a contradiction, by using the linear independency of Lemma [7. II twice 
with respect to the variables u and v for the corresponding curve of (II. ip . Hence, all the 
Aj must be polynomials. Hence, we see that the desired right hand side must be expressed 
by using constants a, b, c, d, e, /, 5^1, g2, hi, /i2, H, i2,j, ki, /c2, which are polynomials in /XjS 
and independent of the Ui and Vi , clS follows: 
(8.3) 

« [Pll{u) - Puiv)] + b [pl2iu)p23iv) - pl2{v)p23iu)] + C [pl3{u) P22{v) - Pl3(f )P22(m)] 
+ d [Ql333(^)p33(f) - <5l333(^0p33(^)] + e/il [pl2 (^^) p33(f ) " Pl2{v)p33{u)] 
+ f[pl3{u)p23{v) - pl3iv)p23iu)] + gi[pl3iu)p33{v) - pig (i;) ^33 (u)] 
+ 92[Ql333iu) - (5l333(^^)] + ^^1 [P23(m) P22 (^^) " P23iv)p22iu)] + h2[pl2iu) - puiv)] 
+ H[p22iu)p33{v) - p22{v)p33{u)\ +22[P13(m) - P\3{v)\ + j [P23(m) p33(^^) " p23{v) p33{u)\ 
+ kx\p22{u) - p22{v)\ + k2\p23{u) - p23{v)\ + A;3[p33(^) " P33{v)\ 

We find by computer using Maple, on substituting the expansion (16. ip up to C13 terms 
of (y{u) into (18.30 . and truncating up to weight 18 in the Ui and fj, that 

a = 6 = c=-l, d = \, e = -i/ii, f = -jj,^, g^ = lu^ ^^^^ 

(8.4) 3 3 3 

g2 = hi = h2 = ii=i2= i = ki = k2 = 0, k3 = g/ig. 

as asserted. In the Maple calculation, it is not necessary to assume the polynomial nature 
of the coefficients as functions of the fij. □ 

Remark 8.2. By applying 

(8 5) lA/^A + A^io 

2 dui \ duj dvj J '^^ 

to l8.H we have —pij{u + v) + pij{u) from the left hand side, and have a rational expression 
of several pij...e{u)s and pij...i{v)s on the right hand side. Hence, we have an algebraic 
addition formulae for pij{u)s. 
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Remark 8.3. By putting v = u — {6,0, 0) and letting 5 — >• 0, we can get a "double-angle" 

(T-formula 

(8.6) 

= -Pni{u) - P112P23 + Pl2{u)pl23{u) - Pn3{u)p22{u) + Pl3{u)pi22{u) 

- |P133(m)<5i333(m) + lp33iu)-£^Ql333iu) + (pii2(m)p33(m) - puiu)pi3s{u)) 

+ /^l (Pll3(^)p23 - Pl3(^)pl23) - I (/Ltl^ " yU2) (pll3(^)p33(^) " Pl3 (^) Pl33 (^^)) 

- |/W8P133(m). 

In the case of the elliptic curve, the corresponding relation is a{2u) = —p'{u)a^(u), whilst 
the corresponding formula for the hyperelliptic genus two curve is given in j4], p. 129. 

9. The second main addition theorem 

The second main addition result applies only in the purely trigonal case (12.11) . using 
the results of Lemma 16.31 The formula is as follows: 

Theorem 9.1. The a-function associated with (12. ip satisfies the following addition for- 
mula on J X J : 

(9.1) + + lc].M„ + |c» ^ ^ 

where 

R{U,V) = -lpl3iu)d3Qi333{v) - f p23(?^) Pll2(t') " |plll(?^) + Jpl22iu)p^^'^\v) 

- Ip222iu)p^''\v) + ^,d3Qr333{u)p^''\v) + '^P333{u)p^''\v) - lfl3P333iu)p^''\v) 
+ llJ^6pl3{u)p333{v) - \tl'9p23iu)p333iv) " |Atl2p333(M)- 

Proof. Our goal is to express 

2) a{u + v)a{u + [C]v)a{u + [C» 



a(u)'^a(v)^ 



using several p functions. Because ([92D belongs to r( J x J, 0{3{{Q^'^^ x J) U ( J x e'^l)))), 
a similar argument to that at the beginning of the proof of Th. 18.11 shows that it must be 
a finite sum of multi-linear forms of the 27 functions in Lemma 17. Ij namely, of the form 

finite sum 

(9.3) C,X,iu)Y,iv), 

3 

where Xj and Yj are any of the functions appearing in the right hand side of the description 
of r'( J, (9(30[^1))) in Lemma [7.1^ and the Cj are polynomial in /ijS. Moreover, (19.21) has 
the following properties: 

LI. As a function on J x J, its weight is (—5) x 3 = —15; 

L2. It is invariant under u 1— [Qu (resp. v [C]"^); 
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L3. It has a pole of order 3 on (Gl^l x J) U (J x e'^]); 

L4. It is invariant under the exchange u ^ v (by Lemma [6.31) . 

Hence, (19 .3^ has the same properties. Thus, we may consider only the functions in our 
basis of r{J, C(36[^])) that have the following corresponding properties: 

Rl. The weight is greater than or equal to (—5) x 3 = —15; 
R2. They are invariant under m i— > [Qu; 
R3. They have poles of order at most 3 on 
There are 12 such functions and they are listed as follows: 

1, Pi3 (weight = -6), (weight = -3), 

pill (weight = -15), pii2 (weight = -12), pi22 (weight = -9), 
p222 (weight = -6), p333 (weight = -3), p^^^^ (weight = -12), 

(weight = -9), (weight = -6), 

<93Ql333 = -6(pl3p333 - Pl33p33) - 3pi22 (weight = -9), 

and the p'*-'^ are defined in (17. 2p . Here the last equality is given by cross-differentiation 
from 9i(53333 using the first of the relations in Lemma [5. II with /ii = /U.2 = /i4 = 0. Since 
(19.21) is an even function, it must be of the form 

(9.4) + .)<^(^ + [Cl^M^ + IC» ^ ^ p 

where 

R{u,v) = aipi3{u)pi22iv) + a2pi3iu) dsQisssiv) + 03^23 (m)pii2(^^) + a^piiiiu) 

+ a5Pi22iu)p^'^^\v) + a6p222(^^)p'^^'(^) + ^7 53<5l333 (^) p'^^^O + ^8^333 («) p'^^' (^) 

+ bipi3{u)p222{v) + b2p23{u)pi22{v) + hp23{u) dsQnss^v) + b^piuiu) 

+ b5p222{u)p^^^\v) + bQPs33{u)p^^'^\v) + £1^13(^)^333(1;) + C2 p23 («) P222 (« ) 
+ C3Pi22(m) + Cips33{u)p^^^\v) + C5 ^3(51333 (m) + C^i ^23 («) p333 (^^) + '^2P222(m) 
+ eip333(M). 

By substituting (16. ip into (19. 4p . and comparing coefficients of different mononomials in 
Ui, Vj, we can find the constants ai, ■ ■ ■ , ei depending on the /x^s. Again, in this lengthy 
Maple calculation, it is not necessary to assume the coefficients are polynomial in the 
/ij. □ 

Remark 9.2. By applying 

(9 5) - ( + + log 

to (19. ip . we obtain algebraic addition formulae for standard Abelian functions, which 
would be interesting to compare with those of Remark (18. 2p . 
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Remark 9.3. By putting v = —u + (5,0,0) into (19. ip . dividing through by 6 and letting 
5 ^ 0, we can get an unusual "shifted" a-forniula of the form 

(na\ (^{u-[(\u)a{u-[C'^]u) A r , so ^ ..so , 

(9.6) 373^^ = 2^ Ci [^i(M)9i/i(M) - /.(M)9i^i(M)] , 



an 



where the fi and the (yfj are the even and odd derivative components respectively of the 
formula in (19. ip . i.e. as given in the following table 



Ci 



1 

2 
_3 
4 

1 

4 

J_ 

12 

"4/^9 



pM(n) 
p["](n) 



9i 



Pl22[U) 
PU2iu) 
Pl22{u) 

d-sQisssiu) 

P333(m) 
P333(^^) 



Ci 



_ 1 
3 

_ 1 
2 

_ 1 
4 

1 
2 

2/^6 
~|/^12 



Pl3(Mj 
1 

p[i2l(n) 
pP2l(n) 
Pisiu) 
1 



53<5i333(m) 
Plll(^) 
P222(^) 
P333(m) 
P333(m) 
P333(m) 



Remark 9.4. In the general elliptic case, there appears to be no formulae corresponding 
to (19. ip and (19.60 . However for the specialized equianharmonic case, where p satisfies 

(p')' = V-^?3, 

it is straightforward to show that 

a{u + v)a{u + C,v)a{u + C^f ) 



and 



a{{l-C)u)a{{l-e)u) 



-\{p\n) + p\v)\ 



These seem to be just the first of a family of multi-term addition formulae on special 
curves with automorphisms, which will be discussed in more detail elsewhere. 
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Appendix A: The fundamental bi-differential 
We write the polynomial f{x, y) in (11.11) that defines the trigonal curve C as 
(A.l) /(x, y) = y^ + p{x)y^ + q{x)y - r{x) 

with 

p{x) = nix + q{x) = ^2^"^ + fJ-bX + Hs, r{x) = x^ + it3X^ + Hex"^ + figx + Hi2. 

We describe explicitly the fundamental non-normalized bi-differential (Klein's fundamen- 
tal 2-form of the second kind) J2{{x,y), {z,w)) in (ll.lSp of the curve for {x,y), {z,w) in 
C defined by f{x,y) = 0. 

Following the scheme described in [1] and applied to trigonal curves in [15], [9] and the 

present paper, one can realize T,{{x,y), {z,w)) explicitly as 

F((x, y), {z, w))dxdz 



(A.2) 



^{{x,y), {z,w)) 



[x - zffy{x,y)fyj{z,w) 
with the polynomial F{{x,y); {z,w)) given by the formula 
F{{x, y), {z, w)) = {wy + Q{x, z)){wy + Q{z, x)) 



(A.3) 



with 



(A.4) 



+ w \ w 



f{x,y) 



T{x,z) ] + y[y 



f{z,w) 



w 



+ T{z,x) -Fo{x,z) 



Q{X, z) = - fl2)xz + (2 - fl5)x - Il8+ III 

T{x, z) = 3yUi2 + {z + 2x)ilq + x{x + 2 z)ilq + Sil^x'^z + p{z)q{x) + x^z^ -|- 2 x'^z 
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The term Fq{x, y) vanishes at /xi = /i4 = and is given by 



Fo(x, y) = 032(0; + z)x^z'^ + C22X^2;^ + C2i(x + z)xz + cnxz + cio{x + z) + cqo, 

-2/^4 - 2/ii^/i2 + /ii"^ + 2/i3/ii 



C32 



C22 



'1, 



C2I = yU6/^l - 2 /il/i4Ai2 + /i3yU4 - yUs/^l^ + 2/ii^/i4, 

Cii = 2 (3yUi^/i4^ + /i6/^4 + /ig/^l - 2;Ui/i4/i5 - yUi^/ig - /i2/i4^) 

ClO = -/i5Ai4^ + AtlyUl2 + 2 yUl/i4^ - 2 /ilyU4/i8 + ;U3/i4, Cqq = 



/i4 (/i4^ + 2 ;Ui2 - 2 /i4/i8) . 



We also remark that the expression flA.SP generahzes the Kleinian 2-polar previously 
derived in the hyperelliptic case [1]. 



Appendix B: Quadratic three-index relations 

A complete list of the known relations quadratic in three-index pij^, up to weight —15, for 
the "purely trigonal" case is given below. Note that with care we can obtain an expression 
such that the highest power on the r.h.s. is no more than cubic. The number in square 
brackets [ ] is the weight. A fuller list for the general (3,4) case is given at [13] . 



P333^ = P23^ + 4pi3 - 4p33p22 + 4p33^, [-6] 
P233p333 = -p22p23 + 4p23p33^ + 2/i3p33^ - 2pi2, [-7] 

P233^ = 4p33p23^ + 4/i3p33p23 + P22^ " |pl333 + 4/i6p33 + 8p33pl3, [-8] 

2 2 2 2 r 1 

P333P223 = 2p33p23 + I^3p33p23 " 2p22 + 3^1333 + 2p33 P22, [-8J 

P223P233 = 2p23^ + 2p33p23p22 + 2/^9 + 4p23pl3 + 2/i6p23 + 2/i3pl3 + 2/i3p23^ 

+ /^3P22P33, [-9] 

P222p333 = -l^3p23^ - 4/i3pl3 + 4/i3p22p33 - 2p23^ " 8pi3p23 + 6p33p23p22 - ^p33pl2, [-9] 

P223^ = 4p23^p22 + 4pii + 4p22pl3 + 4yU6p22 " 4p23pl2 " 4/i3pi2 

2 2 2 4 2 r 1 

+ /Lt3 p33 -4/i6p33 + 3^33^1333 - 8 pi3p33 + 4yU3p23p22, ["lOJ 

Pl33p333 = -2p22pl3 + Pl2p23 + |p33pl333, ["10] 

P233P222 = 2/i3pi2 - |p33pl333 + 2p33p22^ + ^1^&P33^ + 16p33^Pl3 

- 2/i3^p33^ + 4p23pl2 + /^3P23P22 + 2p23^p22, [-10] 

Pl23p333 = 4 p33p23pl3 + 2/i3p33pi3 - 2p22pl2 " |p23pl333 + 2p33^pi2, [-11] 
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P223p222 = 8p33pi3p23 " |/^3pl333 + 4/^9^33 + 4/i6p33p23 + 4/i3p33pl3 
+ 2/i3p22^ - |p23pl333 " P33p23 + 4p23p22^ ["H] 

P233P133 = 2/i9p33 + 2/i3p33pl3 + |p23pl333 + p22pl2, [-H] 

Pl23p233 = 2p33p23pl2 + 2/i3p33pi2 - 2p33pii - |p22pl333 + 2p33p22pl3 
+ 2/i3p23pl3 + 2pi3^ + 2/i6pl3 + 2p23^Pl3 + 2yUi2, [-12] 

P333pl22 = -2p33p22pl3 " /^3p23pl3 - 2pi3p23^ - 6pi3^ - 2/i6pi3 + |p22pl333 
+ 4p33p23pl2 + 2/i3p33pi2 - 2p33pii + 2/ii2 + /i9p23, [-12] 

P223pl33 = 2/i3p23pl3 + 2p23^Pl3 + 2pi3^ + 2/i6pi3 - /i3p33pl2 

- 2p33p22pl3 + |p22pl333 + 2p33pll + 2/ii2, [-12] 
P222^ = -4yU3p33pi2 + 8^33^11 - 4p22pl333 + 24p33p22pl3 

+ 4yU3p23pl3 - 8pi3^ - 4/i9p23 - 8pi3/Lt6 + 4p22^ + 4/i3^pi3 

+ ^13^23^ - 4/i6p23^ - 4/13^33^22 + 16/16^33^22 " 8/ii2, [-12] 
P223P123 = -2p23pll + 2/i9p22 + 2pi3p23p22 + 2/i3p23pl2 + 2pi2p23^ 

+ |/^3P33P1333 - 2/i3pl3p33^ - 2/i9p33^ + 2/i3p22pl3, [-13] 
P133P222 = 4p23p22pl3 - /^3P23P12 " 2p23^Pl2 - |/i3p33pl333 + 4/i3pl3p33^ 

+ 4/i9p33^ + 2/i3p22pl3 + 2p33p22pl2, [-13] 
P122P233 = -/^9p22 + 4pi3pi2 + 2/i6pl2 - |/i3p33pl333 + 2p33p22pl2 

+ 4/i3p33^Pl3 + 4/i9p33^ + 2/i3p23pl2 + 2p23^Pl2 - /i3p22pl3, [-13] 
P333pll3 = -2pl2^ - |pl3pl333 + 6^33^13^ + 2/i6p33pl3 + 2p33^pii - 2/ii2p33 
- /^9P33P23, [-14] 
P133 = 3P13P1333 + P12 -4p33pi3 +4/ii2p33, [-14] 
P223pl22 = -2pllp22 + 4p23p22pl2 + |pl3pl333 + 2pl2^ " 8^33^13^ 

+ /i3(2p22pl2 + 4p33p23pl3) " |/^3p23pl333 + /^3^p33pl3 - |/^6pl333 

+ 4/X9p33p23 + (8/il2 + /i3/^9)p33, [-14] 
P123P222 = 2p22^Pl3 + 2p23p22pl2 - |pl3pl333 " 2pi2^ + 16p33pi3^ 

+ /^3(2p22pl2 - 2p33p23pl3 + |p23pl333) + (8/i6 - 2/i3^)p33pi3 

- 2/i9P33P23, [-14] 

P333pll2 = -2p23pl3^ + 2/i6p23pl3 + 2p33p23pll - 2/ii2p23 - /^9p23^ 

+ |pl2pl333 - 4/i3pl3^ - 4p33pl3pl2, [-15] 
Pll3p233 = 2p23pl3^ + 2p33p23pll - 2/ii2p23 - |pl2pl333 + 2/i3pi3^ 

+ 4p33pl3pl2 - /i9p33p22 + 2//6p33pl2 + 2/i9pi3, [-15] 
P123P133 = 2p23pl3^ - 2//i2p23 + |pl2pl333 + 2//3pl3^ + 2/^9^13, [-15] 
P122P222 = -|pl2pl333 + 8p33pi3pi2 + 4p22^Pl2 - 2/i3pi3^ + 2/i3p33pii - |/i3p22pl333 

+ 4/i3p33pl3p22 + fJ'3^p23pl3 + (-2/i3^ + 4/i6)p33pl2 + (2/i3/i6 - 8/^9)^13 

+ 6/i9p33p22 - 2/i9p23^ - f^3fJ'9p23 - '^fisfiu- [-15] 
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The above equations describe the Jacobi variety as an algebraic variety, see also [S] where a 
general matrix construction is given. By eliminating odd powers with the aid of identities 
such as 

P333^ P233^ — (^333^233)^ = 0, 

one can obtain equations of the Kummer variety, J/ {u [—l\u). 
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